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Abstract 

We introduce the notion of fractional bidromy which is the combina- 
tion of fractional monodromy and bidromy, two recent generalizations of 
Hamiltonian monodromy. We consider the vibrational spectrum of the 
HOC1 molecule which is used as an illustrative example to show the pres- 
ence of nontrivial fractional bidromy. To our knowledge, this is the first 
example of a molecular system where such a generalized monodromy is 
exhibited. 

The description and the understanding of molecular spectra have been a 
long-standing goal in the field of molecular physics, both from experimental and 
theoretical points of view (for recent reviews, see PQH] and references therein). 
Vibrational dynamics of sufficiently rigid polyatomic molecules can be well re- 
produced up to a large fraction of the dissociation threshold by an effective 
Hamiltonian which is obtained either by a fit of parameters to a set of measured 
or calculated energy levels [3] or by the application of canonical perturbation 
theory to an ab initio potential energy surface [3]. An important class of ef- 
fective Hamiltonians are the classically integrable Hamiltonians, which allow to 
simplify the study of the dynamics of the system by the use of constants of 
the motion. Among them, we can distinguish the simplest one, the Dunham 
expansion, which in the absence of strong resonances describes accurately the 
vibrational dynamics at low energy near a minimum of the potential energy 
surface. This effective Hamiltonian can be written as a polynomial expansion 
in terms of the actions of the normal modes, which can be defined globally 
on the whole phase space. Resonant Hamiltonians, i.e., effective Hamiltonians 
with fundamental frequencies in resonance PQ, describe the dynamics at higher 
energies where the coupling between, at least, two degrees of freedom cannot 
be neglected. Even if Hamiltonians of this kind are integrable, i.e., the number 
of constants of the motion is equal to the number of degrees of freedom, their 
classical dynamics may not be globally described by action-angle variables since 
the latter are in general only locally defined. The question that naturally arises 
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is how to detect this feature both from the classical and quantum points of view. 
Indeed, the fact that the actions are only local has a quantum counterpart in 
the joint spectrum of the quantum Hamiltonian as it prevents the existence of 
global quantum numbers [S]. 

In this context, monodromy, which is the simplest topological obstruction to 
the existence of a global set of action-angle variables, has become a useful tool 
both in classical [5J [7J [8] and quantum or semi-classical mechanics [SJ [Sj . First 
discovered and developed by mathematicians, the phenomenon of monodromy 
has been exhibited in a large variety of physical systems extending from atomic 
and molecular ones [§] to purely classical systems [JJ [TU] . Such systems have 
a standard monodromy which is either characterized by an isolated focus-focus 
singularity in the associated bifurcation diagram for the local case or by a sec- 
ond leaf which is glued to the main leaf through a line of bitori for the nonlocal 
situation [3J. Both types of monodromies appear in Fermi resonant systems 
with a non zero angular momentum Recently, different kinds of general- 
ized monodromy, such as fractional monodromy |12] and bidromy [13, 14J, have 
been defined and their presence shown in model Hamiltonian systems (see be- 
low for a concrete definition of these generalizations). The next step in this 
study is the determination of physical systems having such monodromies. It 
is in this spirit that we revisit the analysis of the vibrational dynamics of the 
HOC1 molecule with zero angular momentum. We show the presence in this 
molecular spectrum of nontrivial fractional bidromy which can be viewed as 
the combination of fractional monodromy and bidromy. We first describe the 
corresponding bifurcation diagram, which presents a line of curled tori and a 
line of bitori. Fractional bidromy is defined through a bipath, i.e., a set of two 
loops, which are allowed to cross both lines of curled tori and bitori. This is a 
specificity of generalized monodromies with respect to standard ones for which 
the associated loop lies in the set of regular values of the bifurcation diagram. 
We determine the quantum monodromy matrix for a bipath such that only one 
of its two components crosses the line of curled tori. Conclusion and prospective 
views are given in the last section. 

Vibrational dynamics of HOCl. Several studies have investigated the vibrational 
dynamics of HOCl both from the experimental and theoretical points of view 
(details can be found in |15) and references therein). In particular, although very 
accurate ab initio calculations have been undertaken, it has been shown that 
the use of an effective Hamiltonian allows an original and precise understand- 
ing of the qualitative features of the dynamics j!5j . This effective Hamiltonian 
includes energy levels of the ground electronic state with an energy up to 98 % 
of the dissociation energy. The classical Hamiltonian H is expressed in terms 
of the normal modes coordinates (qi,pi, Q2, P2, 93,^3)- <7i> Q2 and q% are close 
respectively to the Jacobi coordinates (r,-f,R), where r is the OH bond length, 
R the distance from CI to the center of mass G of OH (R is very close to the 
OC1 bond length) and 7 the OGC1 angle (7 = at linear HOCl geometry). H 
can be written as the sum of two terms H = Ho + Hp where the Dunham 
expansion Hp, and the resonant part Hp are respectively equal to 

i—1 i>j 
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The parameters Wj, iy, j/yfe, ^ijfc;, fc, fej and fcij for H0C1 can be found in [15] . 
Since the functions ii = (pf + gf)/2, I = pi + q\ + (p§ + g|)/2 and -ff are 
constants of the motion, the system is integrable. An additional 3:1 resonance 
between the modes 1 and 2 can also be considered but we neglect it in this paper 
since it renders the effective Hamiltonian non integrable. 

Bifurcation diagram. Before describing the classical bifurcation diagram, we 
have to say some words about the quantum problem. The quantization rules 
for nondegenerate vibrations are |15) 



where v\ and P = 2v 2 + v 3 are respectively the number of quanta in the OH 
stretching degree of freedom and the polyad number. v% and P are the quantum 
numbers associated to the classical constants of the motion I\ and /. The 
constant h is an effective Planck constant which can, at least theoretically, be 
modified at will to increase the density of energy levels. This is necessary since 
the notion of quantum Hamiltonian monodromy is only defined rigorously in 
the semi-classical limit where h — > [5]. The value H = 1 corresponds to the 
physical problem |15| . 

From now on, we only consider two degrees of freedom and we assume that 
there is no excitation in the OH stretching, that is v\ = and I\ = 1/2. The 
study would be similar for other values of 1\ . The bifurcation diagram is defined 
as the image of the energy- momentum map EM: (q2iP2-,q 3 ,p 3 ) — > (/ = I, H = 
E). This set is constructed from the determination of the critical values of 
EM. where the differentials all and dH are not linearly independent. Using the 
Liouville- Arnold theorem under suitable conditions, it can be shown that the 
preimage of a regular (i.e., not critical) value of EM is a two-dimensional torus. 
This also means that EM defines a torus-bundle with base space the regular 
values of the image of EM and with generic fiber a torus. The preimage of 
critical values is a critical fiber which can be of different types: point, circle, 
curled torus, bitorus, etc. The topology of the corresponding critical fibers can 
be determined from singular reduction theory [5]. The bifurcation diagram 
of the HOC1 molecule which has been derived in Ref. [IS] is displayed in Fig. 
[TJ Note that the corresponding effective Hamiltonian is of degree larger than 
4 and that this bifurcation diagram cannot be mapped on the catastroph map 
discussed in [31 [TT] for polyad numbers above the second bifurcation. To better 
visualize the phenomenon of monodromy, we plot polyads up to P — 44. P = 38 
is the highest polyad number with an ab initio energy level taken into account 
in the fit. As can be seen in Fig. [TJ the bifurcation diagram has a line of curled 
tori and a line of bitori which intersect at the point C. It is also composed of 
two leaves, called the upper leaf and the lower one, which overlap in the grey 
region of the bifurcation diagram. The two leaves are glued together along the 
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Figure 1: (Color online) Bifurcation diagram of the HOC1 molecule for 1% = 1/2 
as a function of P [see Eq. ((TJ]. The energies of the period orbit [B] characterized 
by P3 = <?3 = have been subtracted from each energy in order to clarify the 
plot (see Ref. [TS] for details). The critical points of £A4 are represented by 
solid and dashed lines. Points of the solid lines lift to points or circles in the 
original phase space, while points of the dashed lines lift to curled tori or bitori. 
A representation of these two singular tori is also given in this figure. The red 
dot C is the intersection point of these two lines. The two paths with arrows in 
black correspond to the bipath used to define the bidromy (see the text). The 
grey region is the zone where the two leaves of the bifurcation diagram overlap. 
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singular line of bitori and along the part of the line of curled tori that lies at 
the right of the point C. The other part of the line of curled tori belongs to the 
upper leaf (see also Fig. [5]which displays the quantum version of this bifurcation 
diagram) . 

The vibrational energies can be obtained by a direct quantum computation 
in each polyad, that is as a function of P. However, this calculation is not 
sufficient to construct the quantum version of the bifurcation diagram and a 
semi-classical analysis is needed to establish the nature of the classical dynamics 
associated to each quantum energy level. For that purpose, we introduce the 
canonical conjugate coordinates (Ik,cpk) (k £ {1,2,3}) which are defined by 
the relations = y/2Ik cos (pk and pk = ~y2Lk sin Note that the polar 
coordinates {Ik,<pk) are not defined if pk = qk = 0. The Hamiltonian H can 
then be expressed in the set of coordinates (J, 9) and (J, ip) where / = 2I2 + 13, 
J = 2I2, = </>3 and ip = 02/2 — 03 with the constraints I > J and J > 0. 
The Hamiltonian H is a function of only /, J and ip. One of the actions of the 
system is / which is global and the other one is given by J — J J 'dip / '{2tt) where 
the integral is calculated along the projection 7 of the flow of H on the space 
(J, ip). J depends on the values of / and E. The regular Bohr-Sommerfeld rules 
state that the semi-classical energies are those which satisfy J = h(n + 1/2) and 
/ = h(P + 3/2). Knowing the leaf to which the loop 7 belongs, we associate 
the same leaf to the corresponding semi-classical energy level. The accuracy 
of the semi-classical energy levels allows us to do the same for each quantum 
energy level and to construct the quantum bifurcation diagram. This diagram 
is displayed in Fig. [2] where we observe in the overlap region of the two leaves 
a superposition of two lattices of points. 

Fractional bidromy. Classical monodromy is the simplest topological obstruction 
to the existence of global action-angle coordinates [6j [7] . Let us consider the 
torus-bundle over the regular values of EM.. Due to the presence of certain 
isolated singular tori such as pinched tori, regular tori are forced to fit together 
with a twist which prevents extending the action-angle variables to the whole 
bundle. The system has then a non trivial monodromy. From a quantum point 
of view, we can also analyze the joint spectrum of the corresponding quantum 
operators and look for the manifestation of classical monodromy in this spectrum 
[5]. The bifurcation diagram becomes a 2-dimensional lattice of points labeled 
by the values of the quantum numbers, the energy E and the polyad number P 
for HOC1. Locally, around a regular value, the lattice is regular in the sense that 
we can find a map which sends this lattice to Z 2 . In order to check the global 
regularity of the spectrum, the method consists in taking a cell, transporting 
continuously this cell along a closed loop and comparing the final cell with the 
initial one. If the two cells are different then the system has quantum (or at 
least semi-classical) monodromy [S]. 

The theory of Hamiltonian monodromy has known recent important de- 
velopments, which resulted in the concepts of fractional monodromy [12] and 
bidromy [T31 [2]. These generalizations are associated to particular singular 
tori, i.e., curled tori for fractional monodromy and bitori for bidromy (see Fig. 
[T]for a representation of these singular tori), and to loops in the bifurcation di- 
agram which are allowed to cross these lines of singularities. In [T2], fractional 
monodromy characterizes a line of critical values corresponding to curled tori 
and ended by a point whose preimage is a pinched torus. If we consider loops 
crossing this line then it can be shown that the notion of Hamiltonian mon- 
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odromy can still be denned in a restrictive way where the monodromy matrix 
has fractional coefficients. The bidromy phenomenon was first introduced in a 
three-degree of freedom system similar to the CO2 molecule [IS] and exhibited 
recently in a general class of two-degree of freedom Hamiltonian systems with 
a bifurcation diagram having a swallowtail structure |14) very close to the one 
encountered in Fig. [TJ except that there is in addition a line of curled tori in Fig. 
[TJ The bidromy matrix is defined through a bipath, i.e., a set of two loops fol- 
lowing each a different leaf of the bifurcation diagram. A bipath appears as the 
only way to generalize the notion of monodromy when the bifurcation diagram 
has a swallowtail structure and a line of bitori. Note the difference between this 
case and the standard non local monodromy [2] which is also characterized 
by a second leaf in the bifurcation diagram but where loops surround the line 
of bitori without crossing it. As for fractional monodromy, a method to cross 
unambiguously the line of bitori can be defined. In the quantum version of these 
generalizations, the size of the fundamental cell has to be increased as discussed 
below in one of the two directions in order to parallel transport the cell across 
the line of critical values. 

Since the system considered in this paper has both a swallowtail structure 
and a line of curled tori, we introduce the notion of fractional bidromy which 
can be viewed as the combination of fractional monodromy and bidromy. As 
displayed in Fig. [TJ we consider a bipath such that one of its two components 
crosses the line of curled tori. In this example, we restrict the determination 
of this generalized monodromy to the quantum case. Let u\ (resp. v\) and 1*2 
(resp. V2) be the two vectors defining the initial (resp. final) cell. The quantum 
monodromy matrix M is the matrix such that 



The vectors U2 and V2 are vertical vectors oriented from the top to the bottom, 
while the vectors u\ and Vi are oriented from left to right. To cross the different 
lines of singularities, the size of the cell has to be increased in the horizontal 
direction for the line of curled tori [T3] and in the vertical direction for the line 
of bitori [T3] . The parallel transport of the cell along the bipath is represented 
in Fig. [5] The cell is broken into two cells when crossing the line of bitori. The 
way this line is crossed is not a priori obvious. This parallel transport can be 
rigorously computed by considering the classical monodromy and its relation 
with quantum monodromy [Tl]. After crossing this line of critical values, the 
two cells are transported along the two parts of the bipath. We finally merge the 
two final cells by adding their basis vectors, which gives the red cell of Fig. [5] 
An analysis of Fig. fallows to deduce the following relations between the initial 
and final vectors: v\ — 2u\ + U2/2 and ^2 = 1*2, which leads to the following 
monodromy matrix M: 



Note that the initial vector u[ — u\ — U2/2 is transformed into 2u[ after one 
loop, which leads in this case to a diagonal monodromy matrix. If we consider 
now a bipath crossing the line of bitori at the right of the point C then the same 
monodromy matrix is obtained although each component of the bipath crosses 
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Figure 2: (Color online) Parallel transport of the cell through the vibrational 
spectrum of HO CI. The energy levels of the upper and lower leaves are respec- 
tively depicted in blue (full) and red (open) dots. The effective value of h is 
taken to be 0.5 for the sake of a clearer illustration. The initial cell is broken 
into two cells when crossing the line of bitori. One of the two cells, the grey 
one, belongs to the lower leaf, while the other belongs to the upper leaf. The 
final cell, the red one, is defined by the addition of the two final cells of the two 
paths. 
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two times the line of curled tori. This shows that the notion of fractional 
bidromy is well defined as it does not depend on where the bipath crosses the 
line of bitori. 

Conclusion. This article proposes, to our knowledge, the first example of a 
molecular system where a generalized monodromy is exhibited. We hope that 
this example will motivate systematic investigations of generalized monodromy 
in the rovibronic spectra of resonant molecular systems. Other two degree of 
freedom molecular systems such as HOBr [TB] are expected to present nontrivial 
generalized monodromies, but this situation is not general. As an example, the 
bidromy phenomenon does not exist in the HCP molecule [3] since in this case 
one cannot define a closed bipath turning around the bifurcation points of the 
bifurcation diagram. 
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